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Recent measurements of the scattering of He and Ne atoms at Rh(110) suggest that these two rare-gas atoms 
measure a qualitatively different surface corrugation: While Ne atom scattering seemingly reflects the electron- 
density undulation of the substrate surface, the scattering potential of He atoms appears to be anticorrugated. An 
understanding of this perplexing result is lacking. In this paper we present density functional theory calculations 
of the interaction potentials of He and Ne with Rh(110). We find that, and explain why, the nature of the 
interaction of the two probe particles is qualitatively different, which implies that the topographies of their 
scattering potentials are indeed anticorrugated. 



PACS numbers:79.20.Rf, 68.35.Bs, 73.20.-r 

The presumed simplicity of the interaction of low en- 
ergy rare-gas atoms with surfaces has strongly promoted 
the use of He atoms as ideal probe particles in scatter- 
ing experiments in order to determine the surface atomic 
geometry and the corrugation of the surface electron- 
density [0J^,§]. Reconstructed surfaces have been suc- 
cessfully investigated and even light adsorbates like hy- 
drogen, which are scarcely seen with other techniques, 
can be located. As a consequence, helium atom scat- 
tering (HAS) is one of the leading tools for the analy- 
sis of structural and vibrational properties of surfaces. 
A widely used form for the interaction potential has 
been derived by Esbjerg and N0rskov B who argue that 
the interaction energy of a He atom and a surface is 
simply proportional to the unperturbed electron density 
of the substrate at the position of the He atom. Al- 
though widely accepted, the credibility of this approach 
has been questioned: For Ni(110) Rieder and Garcia ||] 
reported a serious discrepancy in the corrugation am- 
plitude when comparing their measurements with the 
results of the Esbjerg-N0rskov approach taken together 
with good quality calculations of the surface electron den- 
sity. Annett and Haydock B tried to reconcile the dis- 
agreement by introducing an additional term in the inter- 
action potential. This addition, named anticorrugating 
term, arises from the hybridization between the occupied 
Is orbital of the He atoms and the unoccupied states of 
the metal surface and results in an attractive contribution 
to the potential which should be stronger at on-top po- 
sitions than at bridge sites. Harris and Zaremba j?J crit- 
icized the estimates by Annett and Haydock and argued 
that the anticorrugating term should be by more than an 
order of magnitude smaller than what Annett and Hay- 
dock had evaluated. Harris and Zaremba claimed that 
the discrepancy between theory and experiment lies in an 



improper description of the van der Waals contribution 
and a tendency to over-corrugate the He-surface interac- 
tion potential within the framework of the local-density 
approximation (LDA) of the exchange-correlation inter- 
action. A general consensus on the origin of the effect 
noted by Rieder and Garcia || has not been reached. 

Recent measurements || enforced the interest in this 
important question and in fact raised significant doubts 
about the meaning and interpretation of the important 
HAS method. Rieder and coworkers found unexpected 
differences in the interaction potentials and measured 
corrugations when comparing the scattering of He and 
Ne atoms at surfaces. For Rh (110) and Ni (110) and us- 
ing the Esbjerg-N0rskov approach they concluded that 
the Ne diffraction data reflect the corrugation of the sur- 
face atomic structure and the unperturbed electron den- 
sity. In the case of He atom scattering, however, the 
same type of analysis gave an electron corrugation shifted 
away from the atomic positions: The electron density at 
the short-bridge position appeared in HAS to be higher 
above the surface than the on-top site. Rieder's expla- 
nation, following the arguments of Annett and Haydock 
is that especially at the on-top position the He Is 
orbitals, as well as the Ne 2s orbitals, and the empty 
metal s states hybridize giving rise to an anticorrugating 
contribution. In the case of Ne, however, this contri- 
bution is overcompensated by the repulsive interaction 
between the Ne 2p x ,y orbitals and the metal s states. Se- 
vere doubts about this explanation are in place because it 
assumes that the additional term introduced by Annett 
and Haydock is now even dominating the interaction po- 
tential. 

Obviously, there is profound need for a direct calcula- 
tion of the interaction of a rare-gas atom with a metal 
surface. For this it is important that practically no seri- 
ous constraints on the electronic response of the surface 
and the inertness of the rare-gas atoms are introduced. 
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A related important aspect of such theoretical study is 
that the interaction of a He or Ne atom with a metal sur- 
face is an example of weak physisorption, and a calcula- 
tion represents a critical test of the exchange-correlation 
functional used in ab initio calculations. 

We performed systematic density-functional-theory 
(DFT) calculations exploiting two different function- 
als for the exchange-correlation interaction, namely the 
LDA H and the generalized gradient approximation 
(GGA) If not stated explicitly, the below reported 
results refer to a DFT-GGA calculation. The Kohn- 
Sham equations and energy functionals are evaluated 
self-consistently using the full-potential linear augmented 
plane wave (FP-LAPW) method 00. The Rh(llO) 
surface is treated in terms of a super-cell approach, using 
five layers thick slabs, which are separated by a vacuum 
region of 18 A. The slab thickness is rather small for 
a (110) surface, but because of the weakness and local- 
ization of the interaction it is sufficiently large for the 
present study. The energy cutoff for the LAPW wave 
functions is chosen to be E cut = 15.5 Ry, the muffin 
tin radius -Rmt is 1.24 A, the angular momenta of wave 
functions inside of the muffin tin spheres are taken up 
to Z max = 10. The muffin tin radius for the He and Ne 
atoms is Rmt = 0.9 A. For the potential expansion we 
use a plane- wave cutoff of 70 Ry and a (/, m) representa- 
tion (inside the muffin tin spheres) with Z max = 4. The 
k-integration is performed on an equally spaced mesh of 
88 points in the whole two dimensional surface Brillouin 
zone of a (1 x 1) surface cell. For the evaluation of the 
potential energies of impinging He and Ne atoms we use 
a (1 x 2) surface cell. It is interesting to note that, be- 
cause of the small size of the probe atoms compared to 
Rh and the weakness of the perturbation, a calculation 
with a (1 x 1) cell gives practically the same results. All 
calculations are performed non-relativistically. 

As a first test of the accuracy of the calculations we 
studied the equilibrium structure of Rh bulk and the 
clean Rh(110) surface. The theoretical lattice constant 
(oq 1 = 3.89 A, without accounting for zero point vibra- 
tions) agrees well with that measured at room tempera- 
ture (<2q XP = 3.80 A Jl3| ). We find the first layer to relax 
inwards by Ad^/do = —4.9 %, and the second layer re- 
laxes outwards by Ac^/do = +2.3 %, with do being the 
inter-layer distance in the bulk. These results are in good 
agreement with the values obtained by a LEED analysis 
(Ad 12 /d = -6.8 %, Ad 23 /d Q = +1.9 % 0). 

The interaction between the metal surface and the in- 
coming rare-gas atoms is studied within the adiabatic 
approximation. The substrate geometry is kept frozen in 
the scattering event since the Rh atoms are much heavier. 
We also performed calculations with an adiabatically op- 
timized substrate geometry, which resulted in negligible 
changes: The relaxation of the first layer varied by less 
than 0.1 % if a Ne atom impinges at the on-top position. 

The interaction potential energy was calculated for 
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FIG. 1. Calculated potential energy using DFT-GGA (see 
text) for a He (left) and Ne (right) approaching the on-top 
and short-bridge positions of Rh(110) as a function of the 
distance z from the center of the first surface layer. The insets 
show a magnification of the repulsive part of the potential for 
particle energies used in experiment. Solid lines are guides to 
the eye. 



many positions of the He and Ne atoms (see Fig. [l]). 
At a distance of 6.0 A above the surface we find that 
the atom-substrate interaction is negligible and hence 
use this geometry to define the energy zero of our cal- 
culations. It was also confirmed that at this distances 
it makes no difference where the probe atom is placed 
parallel to the surface. The theoretical results for the in- 
teraction potential energy of He (Fig. |l], left) show clearly 
that the potential energy is anticorrugated with respect 
to the unperturbed substrate electron density and atomic 
structure: At the repulsive part and same distance from 
the surface the energies for the short-bridge geometries 
are higher than those of the on-top geometries. Thus, 
for the whole range of particle energies typically used in 
HAS experiments (~ 20 - 100 meV) the classical turning 
point above the on-top position is closer to the surface 
than that above the short-bridge position. The behavior 
of Ne nearing the Rh(110) surface is qualitatively dif- 
ferent. The potential energy curves for Ne atoms above 
the short-bridge and on-top sites almost coincide, but 
in the repulsive regime there is a clear difference and at 
the same distances from the surface the short-bridge po- 
sition is energetically favored over the on-top position 
(see Fig. [j] right). Thus, according to the calculations 
the corrugation experienced by Ne atom scattering cor- 
responds qualitatively to that of the clean surface elec- 
tron density. This anticorrugation for HAS and "normal" 
corrugation for Ne atom scattering agrees with the exper- 
imental analysis of Rieder, Parschau, and Burg ||. Using 
the data of Fig. |j] we estimate the corrugation amplitude 
along the [110] direction, CuT"; at the particle energies 
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used by Rieder et al. g] as the difference in the classi- 
cal turning point over the on-top and the short-bridge 
positions. We obtain Cio aU ~ —0.06 A for He and 
^scatt ^ +0.04 A for Ne. The comparison of these re- 
sults with those derived in Ref. g , Ci§ att = -0.04 A for 
He and Cfo a " = +0.089 A for Ne, shows that our results 
agree qualitatively and even somehow quantitatively with 
those of the experimental analysis. The quantitative dis- 
agreement may be due to the GGA but it may also be due 
to the fact that the measurements were not performed 
for a clean Rh(110) but a H-covered surface, since the 
adsorbates enabled the identification of the on-top and 
short-bridge positions. Thus, the experimental corru- 
gation amplitude for the clean surface was extrapolated 
from the measured data by assuming that the H atoms 
give rise to a Gaussian contribution to the clean surface 
electron density || . 

Our DFT-GGA results reproduce not only the experi- 
mental corrugation but are also consistent with other fea- 
tures of the probe atoms potential energy. For example, 
the calculated potential well for a He atom is 13 meV 
(8 meV, if we include the zero-point vibration), which 
compares nicely with the value derived from selective ad- 
sorption measurements, 8.2 meV [B. For Ne we find 
18 meV (11 meV with the zero-point vibration). On the 
other hand, it is interesting that with the LDA exchange- 
correlation functional these quantities are in poor accor- 
dance with the experimental data: The turning point for 
He and Ne are systematically closer to the surface and the 
potential wells are too deep (27 meV for He and 61 meV 
for Ne). This is consistent with the well known behavior 
that the LDA systematically gives rise to an overbinding 
in poly-atomic systems. 

In order to analyze the differences of He and Ne atom 
scattering we discuss the changes in the electron den- 
sity of the surface and of the rare-gas atoms induced by 
the interaction. Figure [2] displays the difference between 
the self-consistent electron density of the interacting sys- 
tems and the superposition of the densities of the clean 
Rh(110) surface and a He (Ne) atom. Three positions 
of the probe atoms are selected which correspond to a 
slightly attractive interaction (left panel), to the mini- 
mum of the interaction energy (middle panel) , and to a 
position in the repulsive part close to the turning point 
assuming a kinetic energy of 150 meV (right panel). The 
figure shows clearly that both rare-gas atoms change the 
substrate surface electron density noticeably and that 
they are also significantly polarized themselves. Thus, it 
is obvious that it is not the unperturbed surface electron 
density which is probed by the scattering. Figure || shows 
that for the surface the main changes occur in the d shell. 
The leading effects for both probe atoms are a depletion 
of the d 3z 2_ r 2 states and an increase of d xz and d yz occu- 
pation at the on-top position, and a depletion of the d xy 
states and an increase of d xz and d yz at the short -bridge 
site (The x axis lies along the short-bridge direction and 
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FIG. 2. Density difference plots (see text) for a He (top 
panels) and Ne (bottom panels) atom over the on-top po- 
sition at three different distances from the topmost surface 
layer. The crosses indicate the positions of He (Ne) above the 
surface. The distances are from left to right: 4.2 A, 3.7 A and 
2.2 A for He, and 4.4 A, 3.7 A and 2.4 A for Ne. Full lines 
note an increase and dashed lines a decrease in electron den- 
sity; the values are ±0.5, ±1.0, ±2.0, and 10" 3 bohr" 3 . 

the y direction is parallel to the long-bridge direction of 
the (110) surface). The largest effect happens for the d xz 
electrons: Compared to the unperturbed surface the d xz 
contribution is increased at the turning point by about 
1 %. This increase is partially due to the fact that the 
Pauli repulsion of the rare-gas atoms with the spilling out 
substrate s electrons is reduced by transferring s electron 
density into the d xz band, which for Rh has a particu- 
larly high density of filled and empty states right at the 
Fermi level. We come back to the important role of the 
d xz band below. 

While the reaction of the substrate is similar for both 
probe atoms, the polarization densities of the He and Ne 
atoms are clearly different as is the nature of their inter- 
action with the surface. For He the interaction is medi- 
ated by the He Is electrons and a polarization of the He 
atom away from spherical symmetry which implies a hy- 
bridization of Is and 2p z orbitals, clearly visible in Fig. || 
(top right panel). On the other hand, the interaction of 
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Ne is dominated by the 2p electrons and the easier polar- 
ization of the Ne atom which requires a 2p — > 3s virtual 
transition. This contradicts the interpretation of Rieder 
et al. H which was based on a strong involvement of the 
Ne 2s orbitals. As our calculations show (see also Fig. ||) 
at the on-top position He exhibits a reduction of the Is 
density and an increase of the 2p density (mainly 2p z ); 
a result similar to that found in studies of He physisorp- 
tion at a jellium surface On the other hand, for Ne 
we find a reduction of the 2p z occupation and a slightly 
stronger localization of the p x ,p y states. At the short- 
bridge site the polarization of the He atom is similar to 
that at the on-top geometry, although weaker, but that 
of Ne is different, i.e., here the occupancy of all three 2p 
states is reduced. 

The results are understood as follows. The reflection of 
He and Ne atoms happens rather close to the surface, at 
a distance slightly less than 3 A, i.e., closer than the po- 
sition of the physisorption well. Here the DFT-GGA ap- 
proach describes the interaction with sufficient accuracy. 
The nature of the interaction is determined by electron 
polarizations and hybridizations. Thus, it is not the to- 
tal electron density of the substrate surface which deter- 
mines the interaction, but the electronic wave functions 
which lie close to the Fermi level. The He Is orbital and 
the Ne 2p orbitals interact with these substrate states in 
a qualitatively different manner. For Rh (110) the states 
which are most important at the distance of reflection 
have d xz character, because these substrate states give 
rise to a very flat band (thus high density) which crosses 
the Fermi level close to the Brillouin zone boundary. In 
other words, these orbitals change phase when going from 
one Rh atom to the next one along the short-bridge di- 
rection. Thus, the Bloch state is bonding in character, al- 
though, due to the weak overlap (reflected by the band's 
flatness) one may call it non bonding. 

At the short-bridge position these d xz substrate states 
are thus symmetric with respect to mirror planes along 
xz and yz, and their electron density is low. A He atom, 
with its Is state, feels a Pauli repulsion with these states 
and is efficiently reflected. On the other hand, for Ne 
the 2p x and 2p y states are antisymmetric. Thus, at the 
short-bridge position they will not interact with the d xz 
band. Only the Ne 2p z orbital has the same symmetry, 
but due to its narrow lobe and the low density of the 
substrate states at the short-bridge position, it can still 
approach the surface rather close before the Pauli repul- 
sion becomes important. 

The opposite situation occurs for the on-top geome- 
try were the He Is state is orthogonal to the substrate 
d xz wave functions and can thus approach the surface 
quite close up to the point were the repulsion with the 
energetically lower lying d 3z 2_ r 2 states dominates. Some 
of this repulsion is removed by transferring d 3z 2_ r 2 elec- 
trons into the d xz and d yz bands. For Ne the repulsion at 
the on-top site is much stronger: The 2p z electrons inter- 



act repulsively with the substrate d 3z 2_ r 2 electrons, and 
the Ne 2p x interact repulsively with the d xz electrons. 
The Ne 2p y orbital is found to be affected only little. 

These results imply that the interaction between rare 
gas atoms and a surface is significantly more complicated 
that hitherto assumed: It is not the total electron den- 
sity of the surface which is probed, but the interaction 
is determined by the substrate surface wave functions at 
the Fermi level. Although this is more complicated it 
is also more interesting, because these states are impor- 
tant also for the chemical reactivity of the surface. Our 
explanation of the interaction mechanism has interesting 
consequences. For example, we expect similar "anticor- 
rugation" effects for the d metals which belong to the 
same or a direct neighbor column of the periodic table, 
but for systems with a different band structure we expect 
different effects. 
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